Does the Multiple-scattering Series in the 
Pion-deuteron Scattering Actually Converge? 



A. E. Kudryavtsev 1 , A. I. Romanov 1 ' 2 , V. A. Gani 1 ' 3 

^ State Scientific Center Institute for Theoretical and Experimental Physics, 

Moscow 117218, Russia. 
Department of Physics, National Research Nuclear University MEPhI, 

Moscow 115409, Russia. 
Department of Mathematics, National Research Nuclear University MEPhI, 

Moscow 115409, Russia. 



We demonstrate that the well-known answer for the multiple-scattering series 
(MSS) for a light particle interacting to a pair of static nucleons, calculated in the 
Fixed Centres Approximation (FCA), works well for a wide region of the two-body 
complex scattering length a. However, this approach is not applicable in a narrow 
region surrounding the real positive a half-axis, where the MSS does not converge. 
Simultaneously, for real positive a's the 3-body system forms an infinite set of bound 
states. 

1 Introduction 

The multiple-scattering series (MSS) has played a prominent role in the study of meson- 
nucleus interactions. The sum of all rescattering terms where meson scatters back and 
forth between a pair of static nucleons i.e. in the Fixed Centres Approximation (FCA) 
was derived by L. Foldy in 1945 PQ. It was applied to the ^-scattering by Bruckner 
in 1953 [2]. According to Bruckner, the 7T(i-scattering amplitude is the FCA amplitude 
averaged over the deuteron wave function: 



Here f\(f2) are the 7rN scattering amplitudes on the first (second) nucleon of the deuteron 
and i/idi?) is the deuteron wave function, J \ijj d (r)\ 2 df = 1. Expression ([1]) was first ob- 
tained in ref. [3] by summing the set of the multiple-scattering diagrams. More recently 
the terms of the MSS have been discussed in the context of the effective-field-theory (EFT) 
treatment of 7r-nucleus scattering. The first EFT calculation of the nd scattering ampli- 
tude was performed by Weinberg [JJ. Different aspects of the MSS in EFT formalisms 
have been refined in the next twenty years, see, e.g., refs. [5]. Expression ([1]) for F n d in an 
EFT framework was first recreated in ref. [6]. In our recent publication [7] we discussed 
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from the point of view of an EFT the fact that, if the meson- nucleon scattering is ap- 
proximated by the scattering length, the individual terms of the series are divergent and 
enhanced with respect to the straightforward expectation from the chiral perturbation 
theory (%PT). We also showed in that work that the divergences cancel if the series is 
resummed. 

In this paper we make an attempt to analyse the status of the resummed expression 
for the nd scattering amplitude ([!]). Considering the case of equal amplitudes fi = f 2 = a, 
we get the following expression for the resummed MSS amplitude F(a) in the coordinate 
space at zero energy: 

F(a) = F^(a)+F ( - 2 \a), 
where F^(a) = 2a is the single scattering contribution and 

F^(a) = 2a 2 f |V> d (r)| 2 — *— dr. (2) 
J r - a 

Here F^\a) is the sum of the double-scattering term /^(a), the triple-scattering term 
/®(a) and so on: 

F&(a) = /«(a) + /W(o) + ... + f {n \a) + 

where 

/W(a) = 2a" I \M^\ 2 dr 



Note that we define the scattering length as a = f(k = 0), so the case of a > corresponds 
to effective attraction in the 7riV-subsystem. 

In our previous paper [7J we concentrated mainly on the case of the negative scattering 
length, a < 0. In this case integral (j2]) for the amplitude F^ 2 \a) is well defined. However, 
in the case of attraction (a > 0) the expression for F^ (a) (j2J) becomes formally divergent 
and one needs additional information as to how to deal with the pole singularity in eq. 
(J2]). The first step towards solving this problem has been done in (7J. We showed there 
that the Fourier transform P(Q) of the function (r — a) -1 that enters the integral in (|2j), 

where 

+oo 

/sin x { 7r \ 

■^-^ dx = Ci (y) sin y + ( g ~~ Si ^) cos y ^ 



with Si(y) and Ci(y) being the sine and cosine integral functions, acquires a cut from 
aQ = to aQ = +00, due to the cosine integral having a branch point at y — 0. It means 
that at a > the function P(Q) becomes complex, i.e. it has a nonzero imaginary part. 
Returning to the integral representation for F( 2 \a) (|2J), we conclude that F^(a) has to 
be considered as an analytic function of the complex parameter a = a% + ia 2 - This leads 
us to define the amplitude F^(a) for real and positive a > as the limit of F^ 2 '(ai + ia 2 ) 
as a 2 — > 0. This prescription gives for F^(a) at a > 0: 

F^(a) = 2a 2 V. p. / \^ d (r)\ 2 ^— + 47rW^( a ). (5) 
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According to fl^D , the ird scattering length is complex even when no inelastic channels 
are present. It might mean that the 3-body meson-deuteron scattering problem is not 
well defined if the meson-nucleon off-shell amplitude is approximated by the constant 
scattering length. To check this hypothesis, we can try to find a representation of the 
sum of the MSS terms different from that given by eq. (J2]). 

Our paper is structured as follows. In Section 2 we calculate and resum the MSS 
in the momentum space. In contrast to what we had in the coordinate space, here we 
get a reasonable answer for the amplitude m F^ 2 '(a). However, this answer still reflects 
the presence of the singularity in the expression for F^ 2 '(a) ([2]). The resulting solution 
in the momentum space, obtained numerically, appears to be an oscillating function of 
the size of the discretization grid. To finally decode the answer, we need to introduce a 
finite range into the elementary 7riV-amplitude. This procedure allows us to sum all the 
MSS terms and to get finite answers in both the coordinate and the momentum space 
representations and to perform a comparison of the results. This is done in Section 3. 
In Section 4 we discuss the applicability of both representations of (a) to the cases 
of complex and real negative scattering lengths a. In the Conclusions we summarise 
the results of our research. In particular, we discuss the physical reasons behind the final 
answer for the MSS amplitude becoming ill-defined when the range of the elementary irN 
amplitude goes to zero. 

2 The basic equation for the MSS in the momentum 
space 

Using the Feynman diagram technique one can easily obtain (see, e.g., [3], [8]) 3-dimensional 
expressions for the multiple-scattering amplitudes (a) in the FCA: 



where tpaip) is the deuteron wave function in the momentum space, normalized such that 
/ \ip d (p)\ 2 dp= (2tt) 3 , and 



the IR cutoff x everywhere in the denominators of the pionic propagators. Taking into 
account only the leading S-wave part of the deuteron wave function, one can perform all 
the angular integrations in ([6]), resulting in 




n = 2,3,... 



(6) 
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+oo 

f {2 \a) = 2a 2 ^ J pip d (p){ln K (p,p')}(p d (p')p'dpdp', 
o 

+00 / +00 \ 

J (3) ( a ) = J pMp) I / 2^h*(p,s)ln*(s,j/) <Pd(p')p'dpdp', 



and so on, where we denote 



ln x (j?,p') = In 



(p — p') 2 + x 2 



The answer for the full MSS amplitude m F^(a) = J2 f {n) 

(a) in the momentum space 

n=2 

representation then reads: 

+00 

m F^ 2 \a) = 2a 2 ^- 3 J p Vd (p) ■ R{p,p') ■ Vd (p')p'dpdp' , (9) 


where the function R(p,p') is a solution of the integral equation: 

+00 
f ds 

R(p,p') = lnx(p, P') + a / 7^ ln *(P' s ) R ( s >P')- 

J 2ir 



(10) 



It is convenient to introduce dimensionless variables x, y, z and a parameter L: 

p = p x, p' = p y, s = p z; x = p L. 

To solve equation ({TO]) numerically we introduce the upper limit of the integration 
(regularization), A = Npo. We take the parameter iV to be an integer. To satisfy the 
boundary conditions R(p, 0) = R(0,p') = 0, the parameters Po,x and A have to satisfy 
Po <C x <C A (or 1 <C L < AT). Now we introduce a dimensionless integer- valued grid of 
size A" with the step 1. The discrete version of our basic equation (fT01 reads: 

N 

R N (i,j) = hiL(i,j) + e^2\n L (i,k)R N (kJ), (11) 

k=l 

where we denote R N (i,j) = RN(xi,yj),ha. L (i, j) = \n L (xi,yj),Xi = i - ±, = j - |, 1 < 
i,j < N and e = The solution R N (i,j) of equation ffTTl) depends on the size A" of 

27T 

the grid, and is a matrix N x N. 

Solving eq. ( ITTj) numerically, we found that Rj^(i,j) nontrivially depends on A^Q 
Some examples of the solution Rn{x, y) at different A" are presented in Figs, [la], [Tbl The 
evolution of R^(x,y = 0.5) as a function of A" is the following. At small A" the exact 



1 Note that special efforts were performed to check stability of the results of our numerical procedure. 
In particular, we varied scale po —> po = Po/n and N N — N ■ n (n = 2, 3, 4) and found that numerical 
results survive. 
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Figure la: The dependence of the function R(x, y) on x at y = 0.5 for some typical 
cutoffs N (here iV is the size of the grid). The solid curve 1 corresponds to iV = 50; 2 
- to AT = 150; 3 - to N = 190; 4 - to TV = 200; 5 - to N = 210; 6 - to N = 220. The 
first critical point N^:j t = 205. The dashed-dotted curve shows the Born perturbative 
term R (x,y) = In ^^2+^2 , a l so &t y — 0.5. The parameters used in the calculation: 
a = 10- 3 MeV _1 , p = lOMeV, L = 2.5. 
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R(x,y=0.5) 




x 



Figure lb: The same as in Fig. la but for larger N. Curve 1 corresponds to N = 290; 2 
-to N = 350; 3 - to iV = 450; 4 - to iV = 510; 5 - to iV = 540; 6 - to N = 560. The 

(2) 

second critical point N^ rit = 525. 

solution behaves similarly to the function R (x,y), which is the leading perturbative 
approximation for R^(x,y) in terms of the parameter e, 

R(x v) -\n {x + V)2 + 12 

The function R (x,y = 0.5) is pictured in Fig. [Ta]by the dashed-dotted line. We see that 
the solid curves 1 and 2 in Fig. [la] behave similarly to Ro(x,y). However, with growing 
N the solution Rn(x, 0.5) deviates sharply from Rq(x, 0.5). At N near N^} it = 205, the 
solution blows up to the positive infinity, then changes its sign to become the negative 
infinity and continues its evolution as N increases (compare curves 4, 5 and 6 of Fig. [Tal) . 
Simultaneously, a,t N > N^\ t solution R N (x,0.5) develops a node at small x. This node 
moves to the right with the growth of N (compare curves 5 and 6 of Fig. [Taj. 

The evolution of Rn(x, 0.5) with iV further increasing is shown in Fig. [lb] The solution 

(2) 

blows up again at iV near the second critical point = 525 and the story repeats. At 
iV > N^ it , the second node of Rn(%, 0.5) emerges, see curves 5 and 6 in Fig. [TBI 

The further observed evolution of Rn(%, 0.5) with iV appears to be periodic or cyclic 
- there is an infinite series of critical points N^ it1 and the solution R(x, 0.5) blows up in 
the vicinity of each N^ it . This behavior of the amplitude versus is demonstrated 
in Fig. [2a] for some selected scattering lengths a. Within the accuracy of calculation the 
distance between two adjacent critical points (N^t ~ ^crit) was found to be independent 
of i. This observation means the dependence of N^} it on i is linear, which is demonstrated 
in Fig. [2b] for different scattering lengths a. 

Note that in the vicinity of any N^! it the amplitude F^\a) behaves similarly to the 
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Figure 2a: The dependence of the scattering amplitude m F^ 2 \a) on the grid size N. 
The calculation is performed with the following parameters: i) a = 0.0035MeV _1 ; ii) 
a = 0.005MeV -1 ; iii) a = O.OlMeV -1 . 
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Figure 2b: The dependence of pole positions N^ rit on i for different values of a. Line 1 
corresponds to a = 0.0035MeV~\ 2 - to a = 0.005MeV _1 and 3 - to a = O.OlMeV" 1 . 

situation when a new bound state is created at threshold: it grows to infinity and then 
changes its sign. We can argue that what we observe in our numerical calculation is an 
evidence that, in the limit N — > oo, the 3-body system we are studying has an infinite 
number of bound states. To prove this hypothesis conclusively, we need to search for the 
source of creation of the series of bound states in the 7r<i-system. For this purpose we are 
going to introduce a finite range r$ of the force (which is equivalent to a formfactor) for 
the elementary 7riV-system and to study the dependence of the amplitude F^(a) on Tq. 
This scenario will be the subject of the next section. 

3 Introducing a formfactor into the elementary ttN- 
amplitude 

The elementary 7riV-interaction potential is, in general, characterized by a certain finite 
range of force and hence the t-matrix is momentum dependent. The impact of the inclu- 
sion of a formfactor in the function (a) in the coordinate space was already discussed 
in ref . [7] . Here we want to introduce formfactor in the momentum space and to compare 
the results of both approaches. 

Consider the full off-shell nN amplitude in the separable form: 

f(p,p';E = 0) = ag(p)g(p f ), 

where we take 

( )~ M 

y/M 2 + p 2 ' 

The limit of the zero-range 7riV-interaction corresponds to M going to infinity. 
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In the coordinate space, the introduction of this formfactor leads to the following 
modification of the expression for F^ 2 '(a) ([2]): 

= / ^ r _ { ^ He _„_ e _ M ,.^ . (12) 

In the limit x — > this expression coincides with that discussed in Ref . [7] . In contrast to 

what we had in the point-like limit (j2J), this expression for (a) is well defined in the 
x x 

region Ma < 1 H . At Ma > 1 H the pole singularity in the denominator emerges 

s M M & J 6 

again. 

In the momentum space, the expression for the amplitude m Fff (a) with the formfactor 
g(p) reads: 

m F$(a) = 2a 2 — (^ M2 _ x2 j J p^ d {p)R M {p,P> d {p')p'dpdp' (13) 
where Rm{p,p') is the solution of the following integral equation (cf. eq. ffTUj) ): 



M 2 / ds 

Rm(p,p') = ^ Xj m(p,p') + a-— / — hx XjM (p, s)R M (s,p' 

M z — X z J Z7T 




(14) 



where 

lnx,Af(p,p') = ln*(p,j/) - ln M (p,p')- 

Comparing the functions F^ (a) and "F^ (a), calculated in both the coordinate and 
momentum space representations, we find that they practically coincide for the values of 
the parameter M < M , where M is determined from the equation: 

A/„a=l + -^. (15) 

(2) 

For M > M the function (a) in the coordinate representation is formally divergent 
similar to the point-like case discussed in the previous section. 

The function m F^\a), calculated in the momentum space according to eq. (|T3l) . is well 
defined also for M > M . Some examples of the solution for m Fff(a) in the momentum 
space versus M are demonstrated in Fig. [3] for some different values of M. For M > 
M these functions depend crucially on M. The amplitude m F^(a), calculated in the 
momentum space representation at a = 5 ■ 10~ 3 MeV -1 as a function of M, is shown in 
Fig. HJ One can see the set of critical points M^ it which appear as spikes. All the singular 
points M = M { 2 it are located to the right from the point M = M flTSJ), M^ it > M . At 
M < M the amplitude m F$(a) is smooth, positive and finite. For M -> M c ( r ^ m F^ } (a) 
goes to +oo, then changes its sign and continues its further evolution with M. A similar 
behavior is seen near the second critical point M crit . The function Rm{x,V = 0.5) is 
nodeless for M being in the interval M e (x, M^j t ), while it develops its first node in 
the interval M e (M^ it , M^j t ). Some typical pictures for the function R M (x,y = 0.5) are 
demonstrated in Fig. |5j The further evolution of Rm{x,V = 0.5) with M also appears to 
be cyclic. 

We conclude that for M < M the 3-body system under our study appears to be 
in an attractive but perturbative regime without forming any bound state. At M = 
M^ it the attraction is getting strong just enough to form the first bound state. The 
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Figure 3: The function (a) versus a for different values of M. Curve 1 corresponds 
to M = 150MeV; 2 - to M = 220MeV; 3 - to M = 250MeV; 4 - to M = 300MeV; 5 - to 
M = 500MeV. 

amplitude (a) for M > M^) it becomes negative, which corresponds to the effective 
repulsion. With the further growth of M this picture repeats periodically. In the limit 
M — > +oo an infinite set of bound states in the Tid-system is formed. Hence we conclude 
that the perturbative approach is not applicable in this extremely unperturbative case of 
positive scattering length a. Looking again at Fig. HJ we also conclude that the function 
m F^\a) behaves akin to tanM, i.e., periodically oscillates from minus to plus infinity. It 
means that this function has no definite limit as M goes to infinity, i.e., in the point-like 
limit. This observation is in line with the fact that the function F M (a) in the coordinate 
representation is not well-defined for M > Mq. On the other hand, treating F^(a) as an 
analytic function of a = a\ + ia2, we come to expression (jSJ). 

We thus obtained two different answers for the function F^ 2 \a). The straightforward 
solution of the integral equation for the function Rm(p,p') at real positive values of scat- 
tering length a tells us that the function m F^(a) = lim m F?)(a) does not exist. The 

M— >+oo 

analyticity hypothesis for the function F^(a), on the other hand, gives us a unique solu- 
tion for F^(a) (see eq. (jSJ)). The situation looks rather ambiguous. To make the correct 
choice and to solve the problem we should study the properties of the function F^> (a) in 
the complex a-plane. It will be the subject of discussion in the next section. 

4 Inclusion of absorption 

In the real 7T(i-system the effect of absorption is not the leading one because both elastic 
7riV-amplitudes (n^p — > n~p and n~n — > n~n) are real. The role of absorption is more 
crucial in the K _ <i-scattering problem because of the presence of open inelastic two-body 
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Figure 4: The function m F^'(a) versus M for a = 0.005MeV~\ M = 222MeV. 

channels K~N — > 7rS(7rA). The FCA is widely used for the K~d scattering problem, see, 
e.g., ref. [9]. For more recent publications see ref. [10J and references therein. 

Note that expression (T5]) for the amplitude F^ 2 \a) is well defined only if the imaginary 
part of the scattering length a = a± + ia,2 differs from zero, a,2 7^ 0. On the other hand, 
we may solve the equation ( TT4"|) for the function R M (p,p') considering the parameter a 
as a complex one. Calculating this function m F^\a) in the momentum space, one can 
compare the results of calculations for Fff(a) in both the momentum and the space 
representations. This comparison is demonstrated in Figs. I6T-TT31 The real part of the 
scattering length was fixed, 

a x = 5 ■ lO^MeV" 1 , 

and we varied the imaginary part starting from very small positive 0,2- The starting 
point is Fig. [6] that reproduces the function m F{^(a) versus M for real and positive 
a = ai — 5 • 10~ 3 MeV -:L . The function m F^\a) in this figure is analogous to what is 
shown in Fig. HI We see quite a regular behavior of m Fl { (a) at small M < M and several 
singular points (or spikes) at M > M . The imaginary part of m F^\a) is equal to zero. 

In Fig. we draw the real and imaginary parts of the function ™F« (a), calculated 
in the momentum space for complex scattering length a = 0.005 + 0.00005? (MeV -1 ). 
Both Re m Fj^ } {a) and Im m F^ (a) are shown in this figure. The behavior of the real 
and imaginary parts of m F^\a) at M near each of the M^ it looks similar to their be- 
havior in the vicinity of a Breit-Wigner resonance. The next step of the evolution of the 
function m F^\a) with increasing Ima = a2 is shown in Fig. [Hb (0,2 = 0.0001 MeV -1 ). 
Although for this value of a 2 the singular points at this figure are still separated, their 
corresponding regions already start to overlap. This tendency is better seen in Fig. [9b 
(a 2 = 0.0003 MeV" 1 ) where the isolated singularities are seriously squeezed or deformed. 
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gure 5: The function Rm(x, y = 0.5) versus x at different typical M: a) M — 150MeV; 
M = 260MeV; c) M = 280MeV; d) M = 330MeV; e) M = 380MeV. 
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Figure 6: The function m F^\a) for scattering length a = 0.005MeV _1 . 

Further evolution of m F^ I (a) with the growth of a 2 is demonstrated in Figs. llOb and [Tib. 
We see in these figures that isolated peaks disappear and both the real and imaginary 
parts of the function m F^ (a) become smooth functions — no singularities any more. This 
tendency is confirmed in Figs. [T2b. fT3b drawn for even larger values of a 2 . We observe 
therefore that if the absorption (Ima) is not extremely small, i.e., a 2 is comparable with 
ai, all the calculated curves behave quite regularly. 

Now it is time to compare these results with what we get in the coordinate representa- 
tion, i.e., with the integral representation given by eq. ( TT2|) . In Fig. [7^ the real and imag- 
inary parts of F M '{a), calculated according to eq. (fl2j) for a = 0.005 + 0.00005? (MeV~ ), 
are shown. Both functions look quite regular and smooth. Note that lm.F M (a) is ex- 
tremely small at small M < M w 220. This is a consequence of the fact that Ima = a 2 

( 2) 

is taken extremely small. The growth of Im F\ 4 (a) for M > Mq looks like a threshold 
phenomenon. The origin of this non-negligible imaginary part is similar to what is given 

in the r.h.s. of eq. ([2D for the point-like case. This contribution may be called the "un- 

(2) (2) 
physical" part of Im m F^'{a). In the limit a 2 — > and M — > 00 the function ImF| f (a) 

coincides with the imaginary part given in eq. pj). Further evolution of F M (a) at dif- 
ferent values of a 2 is demonstrated in Figs. [8^1 - [T3h . As it is seen, the contribution of 
the "physical" absorption is getting more prominent with growing a 2 , which is intuitively 
clear as absorption on each individual scatterer is getting stronger. 

It is instructive to compare Figs. [Th^ and [7b, drawn for the same values of a 2 . At 
extremely small a 2 they look absolutely different. The function "^/(a), calculated in 
the momentum space for small a 2 , is an oscillating function of M with singular points. 
Because of the oscillating character no definite answer for this function in the point-like 
limit M — > +00 exists. However, with the growing a 2 m F^ (a) is getting smooth and for 

(2) (2") 

large a 2 the profiles of the functions F^' (a) and m F^' (a) calculated by different methods, 
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coincide with each other (compare, e.g., the curves in Figs. [T3k and fl3b) . Note that at 
relatively large a 2 

lim F$(a) = F^(a), 

where F <y2 \a) is the point-like limit for scattering amplitude given by eq. ([2]). 

Concluding this section, we would like to stress that we used the momentum-space 
representation to check the validity of the simple coordinate-space representation, given 
by eq. ([2]) . We proved that almost everywhere in the complex plane of parameter a one can 
use eq. (j2J) to calculate the MSS. However, this equation gives incorrect answer in a narrow 
region of complex a. This exceptional region is shown in FigJT4l The border of the dashed 
area is chosen under condition that functions F^\a) and m F^'(a), calculated in the two 
different representations, coincide with accuracy better than 1%. Everywhere outside the 
marked area the FCA given by the simple and useful expression ([2]) is applicable. Note 
that in the case of real negative a < 0, the answers for the amplitudes F^ (a) and m F M (a) 
practically coincide identically as well. 



5 Conclusions 

In this paper we demonstrated that the pion-deuteron scattering amplitude F nd , calculated 
in the Fixed Centres Approximation according to eq. (jSJ), is correct and can be used in a 
rather wide region of the complex scattering length a. Simultaneously, we proved that the 
coordinate representation in the FCA (j2J) is getting incorrect and so cannot be used in a 
narrow region surrounding the real positive a half-axis. This region is marked in Fig. [TU 
For a belonging to this region the interaction between the incident light particle and the 
pair of the fixed centres, while being attractive, becomes anomalously strong. Due to this 
anomalously strong attraction, an infinite set of (irNN) bound states is formed. 

An analogous situation, where an infinite set of bound states emerges as well, occurs 
in the problem of the motion of a particle in a singular attractive potential, e.g., 
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Figure 8: Real (dashed curve) and imaginary (solid curve) parts of the functions F^'(a) 
and m F$(a) for the scattering length a = 0.005 + O.OOOlzMeV -1 . 
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Figure 10: Real (dashed curve) and imaginary (solid curve) parts of the functions F 1 
and m F$(a) for the scattering length a = 0.005 + 0.0005iMeV _1 . 
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Figure 11: Real (dashed curve) and imaginary (solid curve) parts of the functions F i 
and m F${a) for the scattering length a = 0.005 + O.OOliMeV" 1 . 
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Figure 12: Real (dashed curve) and imaginary (solid curve) parts of the functions F 1 
and m F$(a) for the scattering length a = 0.005 + 0.003zMeV _1 . 
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Figure 13: Real (dashed curve) and imaginary (solid curve) parts of the functions 
and m F$(a) for the scattering length a = 0.005 + 0.005zMeV _1 . 
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Figure 14: The region of complex parameter a in which the functions F<? (a) and (a) 
don't coincide with accuracy better than 1%. 



V (r) = 

with /3 > and n > 2, see, e.g., [H], §§18, 35. The discrete part of the spectrum for this 
problem is unbounded from below, i.e., the binding energy of the ground state is infinite, 
whereas the size r of the ground state goes to zero. This is the phenomenon of "the 
falling down to the centre". 

Vitaly Efimov showed in 1970 [12] that the three-body system of identical bosons in the 
limit of large two-body scattering length a displays a specific discrete scale invariance. As 
a consequence, there is an infinite number of 3-body bound states and the ratio of binding 
energies of two successive 3-body bound states is approximately 515. This cyclic behavior 
of the 3-body spectrum in the limit a — > oo has been rederived as the leading order 
(LO) prediction of an EFT, see, e.g., [15] . Note that the Efimov effect was demonstrated 
to survive also in a model consisting of two heavy particles and a light one when the 
light-heavy interaction leads to a zero-energy two-body bound state, see [II] . 

In our case, we found a similar (but not identical) cyclic behavior of the spectrum of 
the system of a light particle interacting with a pair of fixed centres, arising when the 
scattering length a, corresponding to the scattering of the light particle on each of the 
fixed centres, is positive. 

To study this problem, we developed a method alternative to the FCA of summing the 
MSS in the momentum space. In this representation we solved the problem numerically. 
To compare the results of both methods, we introduced a finite range of interaction 
ro = M _1 for each of the scatterers. The point-like limit corresponds to the case M — > oo. 
The solutions of eqs. fflO|) and f fl4|) for positive values of the scattering length a confirmed 
the existence of an infinite set of bound states in the (7riViV)-system. Along with that, 
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we found that no unique solution for the ^-scattering length exists for the case of a > 0. 
The function F M (a), which represents the nd scattering length as a function of a and M, 
is an oscillating function of M and hence gives no well-defined answer in the point-like 
limit M — > oo. We conclude therefore that the ird scattering lenght is not well-defined 
for positive values of a. 

It is also instructive to discuss expression (jSJ) for the scattering amplitude F^ 2 \a). As 
seen from eq. ©, F (2) (a) is complex even for real values of a. At the same time the 
imaginary part of F^\a), calculated in the coordinate space according to eq. fl5J) even 
with extremely small values of the imaginary part of a, is also non negligible, see, e.g., 
Figs. El El Obviously, expressions © and (jSJ) violate unitarity! This violation of unitarity 
is a direct consequence of the singularity of the interaction we are dealing with. In the 
case of a singular attractive potential, due to the phenomenon of "falling down to the 
centre" , there is a sink of particles and hence the unitarity violation. We have a similar 
picture in the fixed centres approximation: the interaction is getting supersingular, and 
the unitarity is violated. 

Still, in the physically interesting cases like the K~d scattering the imaginary part of 
the elementary (K~N in the example at hand) scattering length is comparable with the 
real part. In this case the interaction in the 3-body system is not singular any more. As 
shown in Section 4, formula ([2]) is applicable in this case and both the coordinate and 
momentum space approaches to the 3-body scattering problem give practically the same 
answers for the scattering amplitude F^ 2 \a). 

It is interesting, however, to understand better the physical interpretation of this 
"unphysical" imaginary part. In hadronic atoms the imaginary part of the scattering 
length is usually proportional to the width T n ; of the atomic level. In the absence of a 
phyiscal absorption it is difficult to say what this additional width means. 

In the end it is worth to stress that the problem we studied here is a special case 
of the general 3-body scattering problem, which is usually formulated in terms of the 
Faddeev equations. We took the elementary (7rA r )-amplitudes to be constant and found 
that for positive and real scattering lengths a the 3-body system has an infinite number 
of bound states. The introduction of a formfactor (or a finite range of the interaction) 
does not resolve the problem completely: for values of M larger than some critical value 
Mq, M > M , the system also has bound states. 

To solve this problem, we first need to understand what are the minor modifications 
of the elementary amplitudes that will allow us to avoid the singularity that arises in the 
FCA with constant amplitudes. 
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